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Abstract. We develop a well-posedness theory for second order systems in bounded domains where 
boundary phenomena like glancing and surface waves play an important role. Attempts have previ- 
ously been made to write a second order system consisting of n equations as a larger first order system. 
Unfortunately, the resulting first order system consists, in general, of more than 2n equations which 
leads to many complications, such as side conditions which must be satisfied by the solution of the 
larger first order system. Here we will use the theory of pseudo-differential operators combined with 
mode analysis. There are many desirable properties of this approach; 1) The reduction to first order 
systems of pseudo-differential equations poses no difficulty and always gives a system of 2n equations. 
2) We can localize the problem, i.e., it is only necessary to study the Cauchy problem and halfplane 
problems with constant coefficients. 3) The class of problems we can treat is much larger than pre- 
vious approaches based on "integration by parts". 4) The relation between boundary conditions and 
boundary phenomena becomes transparent. 

1991 Mathematics Subject Classification. 35L20,65M30. 
August 27th, 2009, revised August 12th, 2010. 



Introduction 

The theory for first order hyperbolic systems, which was developed with fluid problems in mind, is by now 
rather well understood. It turned out that energy estimates via 'integration by parts' and characteristics are 
the most important ingrediencies in the theory. 

Second order hyperbolic systems often describe problems where wave propagation is dominant. In bounded 
domains this leads to a large number of boundary phenomena like glancing waves and surface waves. Attempts 
have previously been made to write a second order system consisting of n equations as a larger first order system. 
However, boundary phenomena such as glancing and surface waves correspond to generalized eigenvalues which 
are not handled by the theory for first order systems. Furthermore, the resulting first order system often consists 
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of more than 2n equations which leads to many comphcations. In particular, the first order system must in 
general be augmented by side conditions to guarentee that solutions of the first order system satisfy the original 
second order system. 

In this paper we describe a theory for second order hyperbolic systems based on Laplace and Fourier trans- 
form, with particular emphasis on boundary processes corresponding to generalized eigenvalues. Our theory 
uses pseudo-differential operators combined with mode analysis, and builds upon the theory for first order sys- 
tems developed in [UIl] ■ This approach has many desirable properties: 1) Once a second order system has 
been Laplace and Fourier transformed it can always be written as a system of 2n first order pseudo-differential 
equations. Therefore, the theory of [TJ[2] also applies here. 2) We can localize the problem, i.e., it is only neces- 
sary to study the Cauchy problem and halfplane problems with constant coefficients. 3) The class of problems 
we can treat is much larger than previous approaches based on "integration by parts" . 4) The relation between 
boundary conditions and boundary phenomena becomes transparent. 

The remainder of the paper is organized as follows. In section [1] we state the general problem and provide 
some basic definitions. In section [5] we treat in detail the fundamental problem of a single wave equation in 
a half-plane subject to different types of boundary conditions. In section [3] we first study two wave equations 
coupled through the boundary conditions and then outline a theory for the general case of n second order 
wave equations. This theory proves that all essential difficulties already occur for scalar wave equations coupled 
through the boundary conditions. Numerical experiments are presented in section|4l where we study the different 
classes of boundary phenomena for two wave equations coupled through the boundary conditions. 

1. Initial-Boundary Value Problems for second order hyperbolic systems 

1.1. Well posed problems 

In this paper we want to consider second order systems which are of the form 

utt^ Pfi{D)u + F{x,t), t>o, xen, FeC^in), (1) 
in the lialfspace Q — {xi > 0, — oo < xj < oo, j = 2, . . . , r}. Here 

r 

Po{D) = A,Dl + Y,BjDl (2) 

where 

Ai^Al> 0, = B* > 0, 
are n x n constant matrices, u is a vector valued function with n components and we are using the notation 

X = {xi,. . . ,Xr), D = {Di,...,Dr), Dj^d/dxj, 
Ut — du/dt — Dtu, — DjU. 

Ki t — Q we give initial conditions by 

u{x,{)) = fi{x), ut{x,0) = f2{x). 

We are interested in smooth solutions which belong to L^(il) and satisfy, at the boundary F — {xi — 0, —oo < 
Xj < oo,j — 2, . . . ,r} n linearly independent boundary conditions. 

r 

Lu Cm + ^ Cju.,^ = <?, 5 6 CS°(F). (3) 
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Here Co,Cj are constant nx n matrices, Ci is non-sigular and, without loss of generality, we assume it to be 
normalised 

Assumption 1.1. Ci — I. 

To facilitate the use of Laplace transformation in time, we frequently assume that the initial data are 
homogeneous, i.e., /i = /2 = 0. This is however no restriction, since it is always possible to change variables 
in a problem with general initial data such that the initial data becomes homogeneous in the new variable. 
Since the Cauchy problem is well posed (see section [3?2|) we can extend the definition of the forcing and the 
initial data smoothly to the whole of R''(a;) and determine its solution. Then we subtract this solution from the 
halfplane problem and obtain a new halfplane problem where only the boundary data do not vanish. This is a 
very natural procedure because all the difficulties and many physical phenomena arise at the boundary. 

We now introduce some key definitions that classify the problems according to estimates one can achieve. 

Definition 1.2. Consider ([TJ-Q for F = 0, /i = /2 = 0. The problem is called Strongly Boundary Stable if 
there are constants rjQ > 0, and K > which are independent of g, such that for all ry > ?7o > 0, T > 

£ e-'^'[\\ui-MmT) + \\M-,mHOir)) dt<K l\-'^'\\g{-,t)\\lo^r) dt. (4) 

Definition 1.3. The problem (HJ-Q is called Boundary Stable if there are constants r]Q > 0, K > and a > 0, 
which are independent of g, such that for all 77 > ?7oj T > 0, 

f e'^^'M-Mmv) dt<^ r e-^^%i;t)\\W) (5) 
Jo V Jo 

Here ||u||hp denotes the norm composed of the L^-norm of u and all its derivatives up to order p. Thus ^ 
tells us that we "gain" one derivative while ([5]) says that u is as smooth as the data. The constants a, rjQ are 
very important. If rjo = 0, then we can choose ?7 = ^ for every fixed T > 0. This shows that the solution grows 
at most like with time. If rjQ > 0, then there is bounded exponential growth. This can happen when lower 
order terms are present. 

The boundary estimates allow us also to obtain interior estimates. In section [2.21 we will prove 

Theorem 1.4. Consider (]^-(^ with F — 0. If the problem is Boundary Stable, then we obtain interior 
estimates of the form where OllHO(r) replaced by \\u{-, t)||HO(Q) o.iT'd- ol &?/ a > a + 1, respectively. 

If the problem is not Boundary Stable, then it is illposed. 

Since we can always reduce the data such that only g ^ 0, we could restrict ourselves to this case. However, 
we are interested in differential equations with variable coefficients in general domains. Thus we have also to 
discuss the case that _F 7^ 0. In particular, we have to show that the problem is stable against perturbations by 
lower order (first order) terms of the differential equations. 

Definition 1.5. The problem ([T])-([3]) with /i = /2 — is called Strongly Stable if there exists ryo > 0, T > 0, 
K > and a > 0, which are independent of g and F such that, for all rj > rjQ, 




Clearly, if ^ holds, then the problem is Strongly Boundary Stable. For first order systems the classical theory 
(see [HH]) tells us that also the converse is true: If the problem is Strongly Boundary Stable, then it is Strongly 
Stable. As we will see, after Laplace and Fourier transformation we can write our problem again as a first order 
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system which satisfies all the conditions of the classical theory and therefore the results of that theory are also 
valid for second order systems. In particular, the problem is stable against lower order perturbations both for 
the differential equations and the boundary conditions. (See Appendix of [1]). 

Due to physical phenomena like glancing and surface waves, the problems for second order systems are often 
only Boundary Stable. This leads to 

Definition 1.6. We call the problem (IT|)-(I3|) Stable if it is Boundary Stable and if, for g = 0, there exists r/Q > 0, 
K > and a > which arc independent of F such that, for all rj > rjo, 

^ e-'^'(\\u{-,t)\\HHn) + WM-Mmn) + ^)llHO(r)) dt < ^ e-'^'\\F{-,t)\\Hom dt- (7) 
^ ^ V Jo 

If ([7]) holds, then we can obtain an estimate even when g 0. We split the problem into two; one with g = 
and F and the other with g and F = 0. For the first problem we obtain Q and for the other we use 
Theorem 11.41 

In applications there is often a standard energy estimate, which can be obtained by integration by parts 
provided that g = 0. This estimate can be written as 



H■,0)\\l^^n) + \\M■,0)\\Ho^n)+ I ||F(.,T)||2jo(^)dr 

"'0 



(8) 



where the constant K is independent of F. In this case we need only to show that the problem is Boundary 
Stable. 

Theorem 1.7. The problem is Stable if it is Boundary Stable and, for g = 0, the energy estimate 0) holds. 
One might be tempted to replace the requirement (jT]) by the weaker estimate 



-277* 



(\H;t)\\mn) + \H-MW))dt^§; I e-^^nF{-Mlo(n)dt. (9) 

^ /TV ./n 



K 
rj' 

However, the definition is not stable against lower order perturbations. In section [2.31 we will give an example 
which is algebraically Unstable, i.e., with time the solution loses more and more derivatives. 

Definition 1.8. We call the problem (IT])-(l3]) Unstable if the estimate ([7]) does not hold. 

For first order systems the generalization to variable coefficients (and then to quasilinear equations) uses 
the theory of pseudo-differential operators and requires the construction of a symmetrizer, as described in [1], 
which is smooth in all variables. If the problem is Strongly Boundary Stable, then, as we have mentioned above, 
the same construction can be used for second order systems. If the problem is only Boundary Stable, then we 
have to modify the construction. This can be done but is technically somewhat complicated and the details are 
beyond the scope of this paper. However, we will make the result plausible. 

Since the stability against lower order perturbations is crucial for the generalization of our results to systems 
with variable coefficients in general domains, we shall give a proof in section [1.21 

It is also well known that stability against lower order perturbations allows us to use "localization" to 
decompose an initial boundary value problem on a general compact domain into a finite number of problems 
which are either initial values problems in the whole space, or initial boundary value problems in the half space. 
We illustrate the technique with a simple example in one dimension. 

Consider the initial boundary value problem for the wave equation on the strip 

Mtt = M^x, a:: €[0,1], te[0,oo) (10) 

with initial and boundary conditions 

uix,0) ^ fi{x), utix,0) ^ f2{x), Bou{0,t) ^ go{t), Biu{l,t) ^ gi{t), (11) 
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where Bq and Bi are linear first order differential operators. 

A partition of unity of [0, 1] can be chosen as a set of three C°° functions fi{x), ip2{x), ifsix) where ipi is a 
cutoff function 

ipi{x) = l if a; < 1/4, ipi{x)=0 if a; > 1/2. 

Similarily, 

ip3{x)=0 if a; < 1/2, (^3(0;) = 1 if a; > 3/4, 

and 

ip2{x) = 1 - ^l{x) - (psix). 

We now define the functions 

Ui{x,t) = (pi{x)u{x,t), 1 = 1,2,3. 
Clearly u{x, t) — ui{x, t) + U2{x, t) + u^{x, t) for all x G [0, 1] and 

Uitt = fiUtt = ip^iUxx = Uixx + Li, 

where 

Li = -2ipi^{uix + U2x + Us.x) - (fixxiui +U2+ U3) , 

consist only of lower order terms and has the same support as Ui. Defining also fii{x) — ipi{x)fi{x) and 
f2i = Vi{x)f2{x) we obtain that ui solves the half line problem 

uitt^uixx + Li, a; e [0,00), t>Q, 
ui(a;, 0) = /ii(a;), uit(a;, 0) = /2i(x), i?o"i(0, i) = 3o(0- 

U2 solves the initial value problem on the whole line 

U2tt ^ U2xx + L2, U2{x,0) = fi2{x), U2t{x,0) = f22ix), xe(-00,Oo), t > 0, 

and W3 solves the half line problem 

U3tt = U3xx + L3, a; e (-00,1], t>0, 

U3(x,0) = /i3(x), U3t{x,0) ^ f23{x), BiU3{l,t) ^ gi{t). 

If the three problems for ui, U2 and U3 are well posed, then the original problem pil)) - (im) is well posed. 

To treat variable coefficient problems one can invoke what is known as the "principle of frozen coefficients" 
to replace the problem by one with constant coefficients. Heuristically, one can think that if one localizes the 
problem to very small regions then the coefficients of the equation in each region are nearly constant and the 
behavior of the solution is near to that of an equation with constant (frozen) coefficients. The proof of the 
validity of this requires the use of pseudo-differential theory. Here we claim the validity of this principle for our 
problem but do not go into the details. 

1.2. Stability against lower order perturbations 

If the problem is Strongly Boundary Stable, then it is Strongly Stable and therefore stable against lower order 
perturbations. We shall now prove that the corresponding results hold for well posed problems. 

Theorem 1.9. Consider the problem ([iP~(0i for F = 0, g ^ and change the boundary conditions (0i to 

Lu = g + lu, \l\ bounded. (12) 
// the problem is Boundary Stable, then the same is true for the perturbed problem. 
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Proof. We consider lu as part of the data. Then ([5]) becomes 

r e-'^^ui;t)\\W)dt<^(\l\' re-2"*||^(.,0ll^o(r)dt+ e-'^%i;t)\\W)dt) ■ (13) 

Jo V ^ Jo JQ ' 

By choosing 770 such that -^UP < -r, we obtain an estimate also for the perturbed problem. □ 
Theorem 1.10. Consider the •problem with g — and change the differential equations to 

utt ^ Po{D)u + Pi{D)u + F. (14) 
Here Pi{D) is a first order differential operator with bounded coefficients, i.e., 

IIA(^)"IIho(o) <^ill"llHHa)- 
Assume that our problem is Stable. Then the perturbed problem has the same property. 

Proof. In the same way as in Theorem 11.91 we consider Pi {D)u as part of the forcing and choose 770 and a 
sufficiently large. Then the desired estimate follows. □ 

Remark 1.11. One could argue that there are too many definitions. The reason for the concept of boundary 
stabillity is that it gives the easiest test for wellposedness. If the problem is Strongly Boundary Stable then this 
test provides sufficient conditions for wellposedness. If the problem is only Boundary Stable we can use the test to 
find out what kind of boundary phenomena the problem has. This knowledge is crucial in constructing numerical 
methods. In the Boundary Stable case the test provides only necessary conditions; to obtain wellposedness one 
needs to show that the estimate ([7]) also holds. This is always true if there is an energy estimate for homogeneous 
boundary data. 

Remark 1.12. If we do not want to distinguish between Strongly Stable and Stable, we call the problem well 
posed. 



2. A SINGLE WAVE EQUATION 

2.1. A necessary condition for well posedness 

To derive necessary conditions we consider in this section the halfplane problem for the wave equation 

utt = + Uyy + F{x, y,t), a; > 0, -00 < y < 00, t >0 (15) 
with initial conditions, at t = 0, 

u{x,y,0) = fi{x,y), ut{x,y,0) = f2{x,y) 
and one of four types of boundary conditions ata; = 0, — oo<y<oo: 

1) Ut — aux + buy + g, a,b real, |6| < 1, a > 0. 

2) Ux = ihuy + g, b real, 6^0, |6| < 1. 

ON (16) 

3) Ux = g. 

4) Ux — buy + g, b real, 6 7^ 0. 
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The source F and the data fj,g are compatible smooth functions with compact support. We are only interested 
in solutions with bounded L^-norm and therefore we assume 



/oo poo 
/ \u{x,y,t)\'^ dx dy = \\u\\^ < oo for every fixed t. 
-oo Jo 



(17) 



We start with a test to find a necessary condition such that the problem is well posed. In this chapter and 
throughout the rest of the paper, s = ry + denotes a complex number where 77, ^ G M. 

Lemma 2.1. Let F ^ g — 0. The problem ( [j<5)) - l75)) is not well posed if we can find a nontrivial simple wave 
solution of type 

u^e'^+'^^yipix), \\if{x)\\<oo, Res>0. (18) 
Proof. If we have found such a solution, then 

w„ = e""*+'""V(aa;), a > 0. 

is also a solution for any a > 0. Since Res > 0, we can find solutions which grow arbitrarily fast exponentially. 

□ 

We shall now discuss whether there are such solutions. Introducing (IT51) into the homogeneous differential 
equation (ITSI) and homogeneous boundary conditions (|16l) gives us 

fySxx - (s^ +t^^)'^ = 0, \\ip\\<00. (19) 

pyp is an ordinary differential equation with constant coefficients and boundary conditions 

1) sip{Q) = a<^^(0) + hiu)ip{Q), a, b real, |6| < 1, a > 0. 



2) (p^{0) = -buifiiO), 6 real, 6 7^ 0, |6| < 1. 

3) ¥'^(0)=0. 

4) ipx{0) = biujip{0), b real, 6^0. 



(20) 



The general solution of ([191) of the form 

(^(x) =aie«" + a2e-"", (21) 
where ±k are the solutions of the characteristic equation 



- (s^ + w^) = 0, i.e., K = \/7^ + uj^. 

We fix the argument of ^ by 

-TT < arg(s^ + cj^) < TT, argV + = i arg(s^ + w^). 

From the general theory (also proved in Lemma lA.Sl in the appendix) we know that, there is a constant 6 > 
such that 

Ren > S Res. 

Therefore e if and only if cji =0. Introducing (j21l) into the boundary conditions gives us 

1) s = —an + iujb, a, b real, \b\ < 1, a > 0. 

2) K = a;6, 6 real, b ^ 0, \b\ < 1. 

(22) 

3) K = 0. 

4) K — —iujb, b real, 6 7^ 0. 
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Since, by assumption, a > and ReK > 0, there are no solutions of type (IT51) for the first kind of boundary 
condition. It is important to stress here that chosing the the wrong sign for a in the first type of boundary 
condition results into an ill posed problem and no solution can be computed. The second case in (j22p implies 

^2 + ^2^ ^2^2^ i.e., S^^U^iP-l). 

Thus there are no solutions of type ^TE\\ for |6| < 1, b real. As Rew > there are no solution of type (IT5|) for 
the third boundary condition. Finally, since Re k > and b is real, there are no simple wave solutions of type 
(fT8|) for the fourth kind of boundary condition. 
We have proved 

Theorem 2.2. For the boundary conditions il6]) there are no simple wave solutions of type il8\) other than the 
trivial solution u = 0. 

Since (|19|) and (f22|l define eigenvalue problems, we can phrase the theorem also as 

Theorem 2.3. The eigenvalue problems m9\) and have no eigenvalues with Res > 0. 

We shall now introduce the concept of generalized eigenvalues. For that purpose we write in terms of 
normalized variables. 

Definition 2.4. Let s' — iS^Q, uj' = Wq be a fixed point and consider (|19p . (l^^ for ,s' = i£,'Q + ri', uj' ~ lj'q, rj' > 0. 
(i^Q,a;g) is a generalized eigenvalue for a boundary condition if in the limit ?/ — )■ the boundary condition is 
satisfied. 

We now calculate the generalized eigenvalues. By Lemma IA.7[ there are no generalized eigenvalues for 
boundary conditions of type 1). 

For boundary conditions of type 2), we need to consider 




As Re k' > 0, there will be a generalized eigenvalue for boundary condition 2) if and only if bujQ > and 

Since 

the corresponding eigenfunctions are 

They represent surface waves which decay exponentially in x, i.e. in the normal direction away from the bound- 
ary. They are important phenomena in many applications (e.g. elastic wave equations). 
For boundary conditions 3) and 4) , cjq must satisfy the relation 

Co = ±Vl + fe2wo, k'q = -iuj'ab, 

and where the sign in the first relation is chosen so that £,'of"^'o < (because Ren' > 0.) The corresponding 
eigenfunction is 

U = gi"o(fc2:+y)g±iyr+P"|wo|t /•23'j 
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For boundary condition 4) {b ^ 0,) they are oscillatory in x,y,t. For boundary condition 3) {b = 0), they are 
constant normal to the boundary and they are called glancing waves. They are important physical phenomena 
(e.g. Maxwell's equations). 
We collect the results in 

Theorem 2.5. There are no generalized eigenvalues for boundary conditions of type (1). For boundary condition 
(2), (3) and (4) the generalized eigenvalues are given by 



2) $0 = ± ^0: h'q = buj'^ > 0. 

3) IQ^Kl ^'o^O. 



2.2. Reduction to a first order system of pseudo-differential equations 

The estimates obtained in this and subsequent sections are expressed in Fourier-Laplace transformed space. It 
is clear that all these estimates have their counterpart in physical space (such as the estimates in the definitions 
of section [LT|) . To understand the relation between both types of estimates we refer to chapter 7.4 of [2 and 
chapter 10 of [3]. 

We consider ([T5t - p7|) with homogeneous initial data. We Laplace transform the problem with respect to 
t, Fourier transform it with respect to y, and denote the dual variables by s,u;, respectively. For Res > we 
obtain 

Uxx = {s + CO )ii + F , II — u{x, w, s), F = F{x, uj,s), < x < oo, (24) 
with one of the boundary conditions 



1) su = aiix + ibuu + g^ 

2) Ux = —buiu + g^ 

3) Ux = 5, 

4) Ux — ibujii + g, 



(25) 



and ||m(-, uj, s)|p < c». 

Introducing a new variable by 

Ux = y/WT^v. (26) 
we write the Fourier and Laplace transformed system as a first order system 



Ux = Afu + Fo. (27) 

Here 



u = ( " I , 



1 

/2 



with 



The eigenvalues /x of M are /ii = — k' and fi2 = The boundary conditions at a; = become 

1) s'u = av + ibuj'u + g' , a > 0, |6| < 1, 

2) -buj'u + g, \b\ < 1, 6 7^ 0, 

' , (28) 

3) v = g, 

4) i) — ibuj'u + g' , b real, 6 7^ 0, 
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with g' = 

Remark 2.6. We present in this and the following section the easiest way to obtain the estimates at the 
boundary and in the interior of the domain. To generalize these results to variable coefficients pseudo-differential 
theory is needed. The transformations S and T introduced below (see eqns. pOKSp *) need to be smooth in 
the dual variables. The smoothness condition may fail only at the double root of M. In this case the Kreiss' 
symmetrizer is used to get the estimates as explained in 1 . 

We shall now calculate the solution for the case when F — Q, and estimate it on the boundary. 
The eigenvector of M connected with — k' is given by 



(29) 



The transformation 

is, except for a trivial normalization, unitary and transforms M into upper triangular form, i.e., 

^ ^'^^=[ ^=ITW ^ ^ 

Here S, S~^, d are uniformly bounded and depend smoothly on k' . 
Introducing a new variable by 



V J \v 



= (32) 



transforms (l27l) into 



wj ~ I K )\v 



(33) 

As the solution is in L-^, we have u = and also u = u. Thus the boundary conditions become 



1) (s' + an' — ibuj')u ~ g' , 

2) (k' - buj')u = -g' 

' \- , (34) 

3) K u = —g , 

4) (k + ibuj')u = —g' . 

By p2p. these boundary conditions become singular exactly at the generalized eigenvalues. 

Remark 2.7. From all the boundary conditons the first one is the most benign. By Lemma fA.TI we obtain the 
estimate on the boundary 

\u{0,oj,s)f <j^^\g{u:,s)f. (35) 

In this case we gain a derivative on the boundary and the problem is Strongly Boundary Stable. According to 
the classical theory [Tl[2] the problem is Strongly Stable. Moreover, the principle of localization holds and the 
problem can be generalized to variable coefficients and then to quasilinear equations. It is worth noticing here 
that away from generalized eigenvalues, i.e. when the coefficients on the left hand side of ([M)) are strictly away 
from zero, the estimate ([55]) holds also for boundary conditions 2), 3) and 4) and therefore the problem can be 
treated by the classical theory. 

Because of the previous remark, we only need to study the estimates near the generalized eigenvalues. We 
have 
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Theorem 2.8. The problem m5\) - [Il^ with F = and fi = /2 = has a unique solution in which, 
Fourier- Laplace transformed, is given by 



(36) 



For the different boundary conditions sharp estimates follow. For boundary condition 1) and, "away" from 
generalized eigenvalues, for all other boundary conditions the problem is Strongly Boundary Stable and 



1) \u{Q,w,s)\'< 



ry const. ,„,o 
—5\9\'- 



(37) 



Near generalized eigenvalues the estimates for boundary conditions 2), 3) and 4) are 

2) KO,a;,.)p<H^|gp, 

3) |u(0,a;,s)p < const, -j — — < const 



77(|s|2 +a;2)i/2' 



(38) 



4) KO,..,.)p<^|gp, 



and the problem is Boundary Stable. 

Proof. Clearly, as w = and ii ~ u, ([55)1 is the only solution to l\'3'3\i . We need to consider only a neighbourhood 
of the generalized eigenvalues (i^o,ajo). For the second boundary condition = ±\/l — 6^ 0;^, i.e., 

Since k' 7^ at (j^qjWq), we can use Taylor expansion. A simple perturbation caculation shows that the worst 
estimate occurs for = ui' ~ 0. In this case we have, 



(39) 



l^^ol - buj'a 



Thus we have 



m(0,a;,s) < const. — 5- = const. — ^. 

rj' rj'^ 



A similar perturbation calculation gives 

\k' + ibixj'\ 

which gives, for the last boundary condition 



(40) 
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For the third boundary condition on can do better. Lemma A. 4 with 6 = 0, gives us 

This proves the theorem. □ 



Theorem 2.9. When /i = /2 = and F = 0. The unique solution to our problem, described in Theorem 127 
satisfies the following interior estimates. For boundary condition 1 ) and, away from generalized eigenvalues, 
for boundary conditions 2), 3) and 4) 

2 const. |gP , . 

" < TTS 7- (42) 

Close to generalized eigenvalues, for the corresponding boundary conditions, we have 

2) < l-^l' 



(|s|2+^;2)l/2 

4) ilif < 5f!;:i|9p 

Proof. By Theorem [2781 and Lemma A. 2, the solution satifies 

\\uf<^mu:,s)\\ 

Since, by Lemma A. 5, always Rck > ^477, we obtain from ([36)1 that (|42p is valid. 2) follows because, by Theorem 
12.51 Re k' ~ bujQ and then Re k > const. y^jsP + ui'^. 3) follows from Lemma A. 6 and Lemma A. 4, according 
to 

=^ co^st. < const. ^^^^ = < const.— -J. (44) 

Finally 4) corresponds to 4) of (p8| . This proves the theorem. □ 

2.3. Estimates for homogeneous boundary data 

We consider now the problem (l?7l) . (P5)) with g' = and treat only the cases 2), 3) and 4) where there are 
generalized eigenvalues (see Remark l2.7p . 

For 77 = Res > 0, the eigenvalues of M are distinct and therefore we can transform ^27\i to diagonal form by 
the transformation 



Let u,v he defined by 



Then, (l27l) becomes 



u\ / —K 0\ f u\ 1 f—F 



^Jx VO V V^/ 2k' ^/\IfTlJ^ \ F 



(47) 
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with boundary conditions 

2) {k' -buj')u{0,uj,s) = {k' + buj')v{0,uj,s). 

3) u{0,uj,s) ^ij{0,Lu,s). (48) 

4) {k + ibuj')u{0, uj, s) = {k — ibuj')v{0, lu, s). 

The equations (|T7)) are decoupled and as Rck > 0, Lemma A.l gives for all boundary conditions 
\v{0,uj,s)\' < 



|K'|2ReK|s|2+w2 8|K'|2ReK' (|s|2 +a;2)3/2' 
1 II-FII^I 1 \\FP 



4|K'|2|ReK|2 |s|2+a;2 4|K'|2|ReK'|2 (|s|2+t^2)2 



We use the boundary conditions to estimate u{0, w, s). 

Theorem 12.51 tells us that Re k'q w bijQ > in a neighborhood of the generalized eigenvalue connected with 
boundary condition 2). Therefore the perturbation calculation (j39p gives us 



k' + buj' 



u{0,cv,s)\'= , , , \v{0,u;,s)\' < 
K — bui 

^ const. \\F\\ 



2, , ,.n const. Ili^ll 



(|s|2 +^2)3/2 



(50) 



7?2 (|s|2+w2)l/2- 



The interior estimate in this case follows from Lemma A. 2, Rck = Re k'-\/|s|2 + uj'^ and ([50 

|u(0,a;,s)|2 const. ||F||2 



\\u{;u;,s)r< 



< 



RcK (ReK)2 |s|2 +a;2 

const. ||F||2 



(51) 



772 |s|2 + a;2 



Since the transformation T is bounded, pS)) tells us that the estimates p9|) -([51 ]) are also valid for ii, v. Thus 
the problem is Stable. 

For boundary condition 3), the generalized eigenvalue is k' = and therefore, for Res' > 0, we know only 
that Rck' > (5Res'. However, by Lemma A. 6, we have a strong estimate for |K|ReK and the estimates (|49|l 
become 

\^r. ....^^-o-^- \\Fr ^ const. Ili^lP 



,,-/ M,9 const. Ili^lP 

\\v{;u;,s)f I "2 2 - 

By (|48p . the same estimate holds for u(0, w, s). By Lemma A. 2, we obtain the interior estimate 

,, , const. ||i^||2 

||^(-,c.,s)||2<^^ I ^ (53) 
Since it, v satisfy the same estimates, the problem is Stable. 
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For boundary condition 4), the generalized eigenvalue is Kq = — ii^o^, 6 7^ 0, is purely imaginary. Therefore, 
for Res' > 0, we can only use the estimate Rck' > Sij' . Instead of ([50]) we obtain now 

,,9 k' — iuj'b , , const. , -,,9 



If 



Again, the same estimates hold for m, v. The estimates are sharp. Therefore we do not obtain the desired interior 
estimate and the problem is Unstable. We have proved 

Theorem 2.10. For boundary conditions 2) and 3) our problem is Stable, but not for boundary condition 4)- 

Remark 2.11. The estimates obtained for boundary conditions 2) and 3) tell us that u gains one derivative 
with respect to the forcing in the interior of the domain and half a derivative on the boundary. The problem 
can be localized and generalized to variable coefficients and then to quasilinear equations. On the other hand, 
the estimates for the problem with boundary condition 4) show that not even a fractional derivative is gained 
with respect to the forcing. This, for a second order equation, means that one derivative of the solution is lost 
at every reflection on the boundary. The problem can not be localized. We do not pursue this problem but 
illustrate below this bad type of behavior with a simple example: a first order system with a boundary condition 
equivalent to 4). 

An example: Boundary reflection with loss or gain of differentiability. Consider a system of differential 
equations 

ut = ~Ux, vt ^ Vx, < X < 1, t > 0, (54) 

with boundary conditions 

u{0, t) = v.,{Q, t), v{l, t) = t). (55) 

Then 

w = e^(*'=")uo, t; = e^(*+^)?;o, (56) 

is a solution of ([M)) . 

Introducing into ([55]) gives us 

Mo = Xvo, e^uo = —\e~^UQ. (57) 

Thus we obtain a solution of (l54l) . (1551) if 



e'' = -\\ (58) 

Let 

A = A„ = Trm + A„, n = 1, 2, . . . , 

then ([55)1 becomes 

e^^" = TT^n^ - 27rmA„ - . (59) 

has a solution 

A„ w log Tin. 

Therefore the solution ([55)) grows like 

If the initial data can be expanded into a Fourier series 

u(a;,0) =^e^"MA„), 



then (j60p tells us that the solution loses more and more derivatives with time. 
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Now change the boundary conditions (1551) to 

u,(0,i) = f(0,t), v,{l,t)=u{l,t). 

Then we obtain, instead of (fSS]) . 

— Alio = ''^Oi Xe^VQ — e~^UQ, 

i.e., 

e'' = (61) 

Therefore there is no loss of derivatives. 

Geometrically, the two sets of boundary conditions represents two different situations. In the first case any 
wave loses a derivative when reflected at the boundary. In the second case, it gains a derivative. 

3. Second order systems of hyperbolic equations 

3.1. Two wave equations 

In this section we consider two wave equations coupled through the boundary conditions. 

Ultt = Uixx + Uiyy, U2U = U2xx + U2yy, (62) 

on the halfplane a;>0,— (X)<y<oo, for i > 0, with homogeneous initial conditions 

u.i{x,y,Q) ^ Uit{x,y,0) ^0, i = l,2, i = 0, (63) 

and boundary conditions at a; = 0, 

uix + biU2y = gi, U2x + b2Uiy ^ g2: x = 0. (64) 

Here 6i, 62, are real and gj = gj{y,t), j — 1, 2, are smooth functions which arc compatible with the initial data 
(for example, functions that vanish near t = 0). We want to show that our techniques of section [2] can still be 
used to describe the behavior of the solution. 
Fourier and Laplace transform lead to 

uixx - i^^'^ + s'^)ui ^ 0, U2xx - i^^'^ + s'^)u2 = 0, Res>0. (65) 

Thus we obtain solutions that belong to L^ 

Ml = e^'+^'^^-'^^uio, U2 = e''+''^'«'^''u2o, (66) 

where 

K = \/ up- ^ s'^ , for Rc s > 0. 
We recall here that Rck > when Res > 0. The transformed boundary conditions become 



-KMio + iojbiU2Q = gi, 

iUjb2UiQ — KU2Q = 52- 

A simple calculation shows that (|57| has a unique solution 



(67) 



ngi + iubig2 ^§2 + «'^625i 

"in ^ — — jnr^ "20 = — 5- — J7--7— . (68) 

K^+a;^0i02 K^+a;^0i02 
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Det( . 1 "^^M =^2^ 6^2^162^0. (69) 



if and only if 



There is an eigenvalue to the homogeneous problem (|65|) . (|67|) with gi = §2 = 0, if the homogeneous system 
([67|) has a nontrivial solution. By (|67|). this is the case if 

+ uj%ih2 = + uj'^{bih2 + 1) = 0. (70) 

There are five different situations: 

(1) &162 < — 1. By (|70p. there are eigenvalues s with Res > 0. Therefore our problem is not well posed. 

(2) 6162 — —1. Now s = is a generahzed eigenvalue. For s 7^ the solution becomes 

K.gi + iujhig2 Kg2 + 1^6251 /„-, n 

■"10 = 5 , "20 = ^ ■ (71) 

In general, we can not expect that uio, U20 stay bounded for s — > 0. We need to assume that gj — gjtt 
are the second time derivatives of smooth functions. 

(3) —1 < 6162 < 0. By ((66|) . (|70p . we obtain generalized eigenvalues and eigenfunctions if s = ±i^/l + 6162 |w|, 
K = ^|6i&2| w. The solution ()68p is singular at the eigenvalues and uio, "20 have a first order pole. In 
physical space we obtain surface waves whose amplitudes become large for 6162 — ^ —1. 

(4) 6162 = 0. Now s — ±cj, K — determine the generalized eigenvalues and eigenfunctions. The behavior 
is the same as for the Neuman problem. We obtain glancing waves. 

(5) 6162 > 0. In this case the generalized eigenvalues and eigenfunctions are determined by s = + 6162 uj, 
K = zLi\/bib2 oj. The eigenfunctions are oscillatory and the solution behaves like the solution of the wave 
equation in section [2.31 with boundary condition 4). Thus it is Unstable. 

In summary we can state that the problem is well posed for boundary conditions 2)-4) and Unstable 

for boundary condition 5). Also, there can be numerical diffculties if 6162 + 1 is zero or close to zero. 

When we study the estimates for the problem in the cases 3) and 4) we get completely analogous estimates 
as the ones found for a single wave equation. 

Consider the system 

UUt = Uixx + Uiyy - Fi, U2tt = U2xx + U2yy - F2, (72) 

with boundary conditions 

uix + biU2y = 0, U2x + b2Uiy = 0, a; = 0. (73) 
Fourier and Laplace transform gives us 

Ulxx = {s^ + UJ^)U1 + A, U2xx = (s^ + UJ^)U2 + F2, 

and define 

Uix = y]sp~ra^t)l, U2x = V|s|2 V2. 

Thus, the first order form of the problem becomes 



ui\ ^ 1 — 5- ( 1 \ / ui\ , 1 ( Q 



V|s|2 +UJ-' ' 



(74) 



with boundary conditions 



-\/|s|2 + u!^ V20 + iojb2Uio — 0. 
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We transform (I74[) to diagonal form. Let 



Then 



1 1 



J = 1,2. 



ui\ f ^1"^ \ ui\ 1 ~ 



X / X / X (75) 



with boundary conditions 



-KUlo + jw6i7i20 = -«:wio, , ^ 

... . . 76 

lUjb2UxQ - KU20 = -'KV20- 



In the same way as in section [^751 we can determine —KVio, —KV20 by solving the equations ([75]) for vi,V2 
and reduce the problem ([7 ^ . ([75)) to the previous problem Thus we obtain the same estimates. 

3.2. General systems. The Cauchy problem 

We consider the Cauchy problem for the homogeneous system of the introduction 

utt = Po{D)u, t>0, X e M'' : -00 < < 00, j = 1, 2, . . . r (77) 

were 

r 

Po{D)^A,DI+J2b,DI (78) 

Ai — Al > 0, Bj ^ B* > 0, are positive definite symmetric n x n matrices and u is a vector valued function 
with n components. 

At t = we give initial data 

u{x,0) = h{x), ut{x,0) = f2{x). (79) 
Also, /i, /2 are smooth functions with compact support. 

We want to show that the problem is well posed. Fourier transform with respect to x gives us 

r 

Utt = -Pq{uj)u, Poii^) = AiLul + Bjuj], Po{uj) = P^iuj) > 

i=2 



and 

Introducing a new variable 
gives us 

Therefore we obtain 



u(w,0) = /i(tj), Ui(w,0) = 72(0;), 



51/2- 
ut = Pq' V, 



u\ _ [ p^^^\ fit 



d_ 

di 
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I.e. 

\\u{;t)f + M;t)f = \\ui;0)f + M;0)f- 
This energy estimate shows that the Cauchy problem is weU posed. 

3.3. The resolvent equation 

Consider the Cauchy problem for the inhomogeneous system ([77| with F{x, t) € Cjf^ and with homogeneous 
initial data /i — f2 = 0. 

Uu = Pq{D)u + F. (80) 
Fourier transform with respect to x and Laplace transform with respect to time gives us the resolvent equation 

(^s'^I +\Lj\^Po{Lj')ju^ F, s = i^ + rj, r]>0, u' = uj/\uj\. (81) 

Since Pq — Pg* > 0, there is a unitary transformation which transforms (j81|) to diagonal form 

(^s^ + \uj\'^fij^Uj = Fj, j = 1,2, ... ,n, u — Siu, F = SiF, 

i.e., 

(s + i\uj\^j){s - i\uj\^j)uj = Fj. 
Without restriction we can assume that ^ > 0. Then 



\s + z|w| Vm7| = y'|e+|w|VM7P + 7y2 > ^|s|2 + |u;|2^,, 
\s-i\uj\^j\ = y^lC- + >V = Res. 

Therefore 



\F 



" - V|s|2 + |w|2^^-Res' 

Choosing Ims = i\uj\y^ shows that the estimate is sharp. 
We have proved 

Theorem 3.1. There is a constant K which does not depend on uj such that the resolvent estimate 

|u(w,s)|<^ (82) 

^|s|2 + |w|2Res 

holds. 

The last estimate shows that we 'gain' one derivative, i.e., if the forcing S H^', then the solution € H''+^. 
Therefore we can prove that the Cauchy problem is stable against lower order perturbations. 

Theorem 3.2. Consider the Cauchy problem with homogeneous initial data for 

utt = {Po{D)+Pi{D))u + F. 

Here Pi{D) represents a general first order operator. There is an rjo > such that the estimate i82\) holds for 
rj = Ke s > r/Q. 
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Proof. We consider Pi{D)u as part of the forcing. Then gives us 

K \Pi{iuj,s)u\ K \F\ 



\u{uj,s)\ < 



Res^|s|2 + |^|2 Res^ 



Since ^f^^'"'''''^ is uniformly bounded, we choose 770 such that 

K \Pi(iuj,s)u\ 1,., 

Then the desired estimate follows. □ 

We can write the resolvent equation also as a first order system. We Fourier transform ((80|) with respect to 
X- = {x2, . • . , Xr) and Laplace transform it with respect to t. Then we obtain 

r 

AiDju^ {s^I + B{uj_))u- F, B{uj^) ^Y^BjUj^, (83) 

u = u(xi,a;_, s) £ L^(0 < Xi < 00). 
Since Ai > 0, there is a constant ct > such that Ai = Al > al > 0. Introducing a new variable by 



we obtain the first order system 



D,rUM{s,.^)r\+ J '^], (84) 



where 



VkP + l^-P^' 

M = Af(s,W_) = ( A-^s^I+B{u_)) ^ 

The eigenvalues k of Af are solutions of 

AiK'^ipo^ {B{uj^)+s'^l)ipo. (85) 

Lemma 3.3. For Res > 0, there are no Kj with Re^j = 0. Also, there are exactly n eigenvalues, counted 
according to their multiplicity with Re k < and, therefore, n eigenvalues with Re k > 0. 



Proof. Assume there exists a. k — icji which is purely imaginary. Then, by (|8ip . 

(s2/-Po(«wi,«c^-))<Po = (86) 

has a nontrivial solution, i.e., is an eigenvalue of PQ{iuj). Po{ii^) < implies that is real and negative which 
is a contradiction with Res > 0. 

The solutions of ([55]) are continuous functions of w„ . Therefore the number of k with Re k < does not 
depend on cj_ and we can assume that aj_ = 0. Then (1851) reduces to 



{s^I-AiK^)ipo^O. (87) 
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Since, by assumption, Ai has positive eigenvalues fij and a complete system of eigenvectors, we can transform 
([87)) into n scalar equations 

(s^ - ^jK^)^^ = 0, j = l,2,...,n, 

i.e., 

K — ±s/^/JIj, /ij > 0, j — 1,2, . . .n, Re s > 0. 
This proves the lemma. □ 



By Schur's lemma, there exists a unitary transformation U = [/(s,w_) such that 

U*{s,uj_)M{s,uj.)U{s,u^) = 



Mn Mi2 
M22 

where the eigenvalues Kji, Kj2 of Mn and Af22 satisfy Re Kji < 0, Re Kj2 > 0, respectively, for Re s > 0. Clearly, 
the transformed equation (j84l) can be solved uniquely for Res > 0. 

Using ([5^ . we shall now derive estimates for the solutions of ([M| . To accomplish this we consider a more 
general forcing. We replace 

by f , K- = K(xi,a;_,s). 



We Fourier transform (j84p with respect to xi and consider 



iwiM = v|sp + |cZrp{; - Fi, 

IWiU = U — i'2. 

Eliminating v gives us 

Therefore, by (l82t . we obtain the estimate 



^|s|2 + |w|2Res 

Eliminating w, we obtain the same estimate for v. Therefore we have proved 
Lemma 3.4. There exists a constant K > such that, for all lui,uj^, s, 

2K 



{M{s,uj^)~iujj) ' < (89) 



In particular, the eigenvalues k of M{s,u!-) satisfy 



2K 

Using scaled variables 



Re^|>^. (90) 



(91) 
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we can write in the form 

\{M'{s\u^'_)-zJj)-'\<^, (92) 
|Re^'|>^. (93) 

3.4. Reduction to a first order system of pseudo-differential equations 

Now we consider the general halfplane problem ([T])-® with homogeneous initial data fi{x) = fiix) ~ 0, 
coupled to the boundary conditions ([3]). We Fourier transform the problem with respect to x_ — {x2^ X3, . . . x^) 
and Laplace transform it with respect to t and obtain (j83p coupled to the boundary condition 

r 

Ux + {Cos + '^CjUJj)u^ g, u ^ u{0,uj-,s), g^g[uj-,s) (94) 

As in section[2]we have to assume that there are no simple wave solutions for Re s > i.e., that the eigenvalue 
problem consisting of the homogeneous equations (|83)) and (|94p have no eigenvalues s with Re s > 0, otherwise 
the problem is not well posed. 

Now we introduce new variables by 

u^A;^''^ui, Diu = yjlpTl^Zfw (95) 
and obtain a normalised version of (l84l) 

ii 

V 



H{s', J) = A-^'^ {s'^I + B{J_)) M^(^ 



where 

/ 

,H 



and s' = s/-\/|sP + oj'_ — w_ / y^sp + lcJTp are scaled variables. 

Instead of (|85|) we obtain now 

1p — k' Lp, 

H{s' , uj')(p = Kip, 



(97) 



i.e., 

n'^^^A-"\Pl + B{J,))A-'''^, k'= " (98) 

vlsp + k-r 

([98|) is a normalized form of (j85|) . 

Lemmas 13.31 and p.4p are crucial because they garantee that we can use the classical theory (see lemmas 2.1 
and 2.2 in [1]). In particular, as in section[2l away from any eigenvalue or generalised eigenvalue s, with Re s > 0, 
the strong estimate of definition 11.51 holds . Thus we need only to consider the estimates in the neigbourhood of 
generalised eigenvalues. 

Also, for Re s > we can use Schur's lemma to transform M to the upper triangular form ()88p separating 
the eigenvalues with Re Kj < and Re Kj > respectively. Then we use the technique of section [5] to estimate 
the solutions. This becomes particularly simple if there is a standard energy estimate and we need only to show 
that the problem is Boundary Stable. 
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Finally, we want to show that in the neighbourhood of generalised eigenvalues our problem behaves like wave 
equations. 

We shall now derive a normal form of M for Sq — i^Q. We have 



Lemma 3.5. Let 



i^Q, Lj'_^ he a generalised eigenvalue. iJ(tjQ,i^g) is symmetric and there is a unitary 



transformation U such that 



where 



U*HU 



Vo 



Hii 
H22 



4' > 4' > • ■ • > C > 0, 



(99) 



f2 



H: 
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^m+2 



V 



(100) 



(101) 



where k'^ — ±K^ (a;'_g, z^q) are the eigenvalues of I198\). n'^+i ~ if and only if 



Co = /3j|w-oP a'^ eigenvalue of B{oj'_^q), < /3mi„ < < /^^aa;, j = 1,2, 



(102) 



if < /3mm|wl.oP> ^^en > 0. 

Proof. Since i/(a;'_g, i^p) is symmetric, we obtain ((99|) - (jl01l) . 

If Co^ > /^maiclw^Qp, then H is negative definite and all k'^ < 0. Correspondingly, if ^q'^ < /3min|wLoPi then 
all Kf >0. If = 0, then ((Ml) becomes 

(S(c.Lo) - ^ = 0, 

i.e., ^o'^ must be an eigenvalue of B{uj'_q). Clearly, the reverse is also true. If ^q'^ is an eigenvalue of B{u!'_q), 
then there is a k'^ = 0. This proves the lemma. □ 

To simplify the arguments we shall make a strong assumption which we shall relax at the end of the section. 

Assumption 3.6. The eigenvalues k!J are distinct. 

In this case we can choose U = U{uj'_,0 as a smooth function of uj'_,(,' in a neighborhood of a;'_Q,^Q. Also 
there is a constant do > such that, in the whole neighborhood. 



\k'^ ~ — "^0 for all i,j with i 7^ j. 
Finally, we make the perturbation s' = i^' + rj' , —rj'^ < rj' < t/'q, rjQ <C 1. We have 



(103) 



(104) 



Since U*A is strictly positive definite, its diagonal elements ajj > are positive. By using well known 
algebraic results (Gershgorin's theorem) this gives us 
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Lemma 3.7. For sufficiently small rjg which depends only on A ^ and m03\) . there exists a smooth nonsingular 
transformation S = I + rj' Si{uj'_,s') such that 



22 



'an 



H{u;'_,s') = S-'U*HUS= {^^' J + S^CV | ■-. |, (105) 







where djj = Ojj + 0{ri') > and Hn, H22 ire given as before but with distinct eigenvalues. 

Now we can construct the normal form for the resolvent equation (|96p . We introduce new variables by 

and after a permutation we obtain 

Theorem 3.8. In a neighborhood ofuj'_Q,^Q the resolvent equation I196\} can be transformed smoothly into 



dAI] = vwtj^{^' i](f)+ J Vi (106) 



^H{u;'_,s') Oj \ij ^|s|2 + |a;_|2 
where F — —S^^U*aY^F. By il05\) . the system il06\) is composed o/2 x 2 systems 



1\ fu 



where 

Clearly, the 2x2 blocks have the form ((27)) of section [2] 

There are no difficulties to generalize the results to the case that the k| have constant multiplicity (this was 
done in [5] for first orders systems). 

4. Numerical experiments 

In this section we numerically solve the strip problem for the scalar wave equation 

utt^u^x+uyy + F{x,y,t), 0<x<l, 0<y<l, t>0, (108) 
with 1-periodic solutions in the y-direction, 

u{x, y, t) = u{x, y + l,t), < x < I, t > 0, (109) 

subject to initial conditions, 

u{x,y,0)^ fi{x,y), Ut{x,y,0) ^ f2{x,y), 0<x<l,0<y<l, (110) 
and boundary conditions 

u,-buy^ go{y,t), x^O, 0<y< 1, t>0, (111) 
u,=gi{y,t), a; = l, 0<y<l, <>0. (112) 
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Here 6 is a constant. We are interested in the three cases b = 0, b real, and b purely imaginary, i.e., 

b = i/3, (5 real. 

To solve the latter problem using real arithmetic, we introduce real-valued functions u*^^^ and m^^^ such that 

w = u(i)+m(2). (113) 
Inserting (|113l) into (|108p leads to the system of scalar wave equations, 

=uW+uW+ReF(x,2/,t), 0<x<l, 0<y<l, t>0, (114) 
u[f =u<-^}+u<^J +lmFix,y,t), < x < 1, 0<y<l, t>0. (115) 

Boundary condition piip can be written as 

= Re.go(y,i), x = 0, < y < 1, t > 0, (116) 
=Imgo(2;,i), x^O, 0<y<l, t>0, (117) 



which is of the form with 6i = /3 and 62 = — 
We introduce a grid with grid size h = 1/{N — 1), 

Xj = ij~l)h, .7 ==0,l,...,iV + l, yk = kh, k = 0,l,...,N. 

Time is discretized on a uniform grid with time step 6t > 0, i„ = n6t, n — 0, 1, 2, . . . and we denote a grid 
function by 

'"Ik = vixj,yk,tn). 
The standard divided difference operators are defined by 

D+,vl, = D.^vl, = Do. = I + ^-.) , 

with corresponding notations in the y- and ^-directions. 
Consider the difference approximation 

D+tD^tvlk = [D+.D^, + D+yD^y)vlk + F{xj,yk,tn), (118) 

subject to boundary conditions 

vlo-vl^_^=Q, j-0,l,...,iV + l, (119) 

^In-^Ii=^, j = 0,l,...,iV + l, (120) 

i5o.<,fe-6^oX,fe =.9o(yfe,tn), fc= l,2,...,iV-l, (121) 

Do.vl^k=9l{Vk:tn): fc = 1 , 2, . . . , iV - 1 , (122) 

for n = — 1, 0, 1, . . ., and initial conditions, 

vlk = h{xj,yk), vjl = f^{xj,yk), j = 0, 1, . . . , TV -fl, fc = 0, 1, . . . , TV. (123) 

When b = 0, the difference approximation (|118l) - (|123p satisfies a discrete energy estimate, under the Courant 
time step restriction 

St < Ch, 
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and is therefore stable. The energy method can not be used to show stabihty of the difference approximation 
when b is non-zero and real, or purely imaginary. However, as we shall see below, our practical experience 
indicates that the approximation is stable also when b is purely imaginary. When b is non-zero and real, the 
continuous problem is Unstable. In this case, the difference approximation is convergent for short times. 

To test the accuracy of the numerical solution, we choose the forcing functions F, g^^ and fk such that the 
exact solution becomes the traveling wave 

Uu,{x, y, t) — sin(27r(a; — t)) sin(27ry). 

For the case b = i P, we use the exact solution 

ul^\x, y, t) — sin(27r(a; — t)) sin(27ry), 
ul^\x, y, t) = cos(27r(x — t)) cos(27ry). 

Table [1] shows the max norm of the error Uw — vdX different times for grid sizes h = 10^^, 5 x 10^^, 2.5 x 10"'^. 
All calculations used St — 0.5/i. The error decreases as 0{h'^) for all three values of b, both at time t ~ 1 and 
t = 10. 



Case 


h 


\\Uw - v\\oo{t 


= 1) 


- v\\oo{t 


= 10) 


6 = 


1 X IQ-^ 


7.09 X 10- 


4 


4.69 X 10" 


-4 




5 X 10-3 


1.76 X 10- 


4 


1.18 X 10" 


-4 




2.5 X 10-3 


4.42 X 10- 


5 


2.97 X 10" 


-5 


6 = 0.5 


1 X 10-^ 


7.09 X 10- 


4 


2.95 X 10" 


-'I 




5 X 10^3 


1.76 X 10- 


4 


7.37 X 10- 


-3 




2.5 X 10^3 


4.42 X 10- 


5 


1.84 X 10" 


-3 


6 = iO.5 


1 X 10-^ 


1.44 X 10" 


y 


1.61 X 10- 


-3 




5 X 10-3 


3.61 X 10- 


4 


4.02 X 10- 


-4 




2.5 X 10-3 


9.03 X 10- 


5 


1.00 X 10" 


-4 



Table 1 . Max error in the solution as function of the grid size with a traveling wave exact solution. 

To illustrate how the 6-coefficient in the boundary condition influences the solution, we study the evolution 
from an initial Gaussian pulse, 

with homogeneous interior and boundary forcing functions, F — 0, go — 0, gi — 0. In these calculations we use 
the grid size h = 5 x 10~^ and time step St = 2.5 x 10-3. The evolution between times t = 0.25 and t = 1.25 
for the three cases 6 = 0, 6 = 0.5, and 6 = i 0.5 is shown in Figure [T] The solution initially propagates outwards 
towards the boundary (first column). Before the pulse reaches the boundary, the solutions are identical for 
all three cases. The middle column in Figure [T] shows the solutions after the initial pulse has reached the 
boundary, at t — 0.75. At this time, there are only minor differences between the three solutions. The influence 
of the boundary condition is becoming more obvious in the right column, corresponding to t — 1.25. While 
the differences between 6 = and 6 = iO.5 are still small and located near the left boundary, the case b — 0.5 
has developed a structure near that boundary which is not present for 6 = 0, or 6 = iO.5. For later times, the 
boundary structure develops into a diagonal streak which extends further and further into the domain until it 
gets reflected by the Neumann condition on the right boundary, see Figure [51 At later times, the reflected streak 
develops a new streak which grows into the domain in the same way. The solution is eventually dominated by 
these streaks which appear to propagate in the direction y + 0.5x. Note that the propagation direction y -I- 0.5a; 
is consistent with eigenfunction (|23|) since 6 = 0.5, see case 4) in section [2?T] 
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0.5 1 0.5 1 0.5 1 




0.5 1 0.5 1 0.5 1 



0.5 








0.5 1 0.5 1 0.5 1 



Figure 1. The solution at times t = 0.25 (left column), t = 0.75 (middle column), and t = 1.25 
(right column) for 6 = (top row), b = 0.5 (middle row), and b = iO.b (bottom row). The 
bottom row is showing the real part of the solution (u*^^'). 



It is also interesting to monitor the max norm of the solution for longer times when b = 0.5, see Figure [3l Note 
that the solution grows exponentially with time, illustrating the Unstable nature of this boundary condition. 
Also note that the solution is slightly larger on the finer grid. This behavior agrees with the predicted exponential 
growth proportional to because higher values of are captured on the finer grid. Note, however, that 
this growth is not due to numerical instabilities because the accuracy test shows second order convergence, at 
least up to i = 10, see Table [T] 

To more clearly see the difference between the cases b — and b — i /3 we take F = 0, go = gi — Q and change 
the initial data to trigger a surface wave, 

fi(x,y) = Us{x,y,0), f2{x,y) = Us{x,y,-6t), 

where 

u,(x,2/,t) = e-l'''^°l^ [cos(t^o(y- Vl^O) +*sin(wo(y-yT^i))] , /3wo > 0. (124) 

This wave decays exponentially away from the a; = boundary with a harmonic oscillation in y, see Figure 21 
The surface wave propagates in the positive y-direction with a wave speed proportional to ^/l~^/3^. As /3 — ?► 0, 
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0.5 1 0.5 1 0.5 1 




0.5 1 0.5 1 0.5 1 




Figure 2. The solution of the unstable case {b — 0.5) at times t = 2.5 to t = 4.5 in increments 
of 0.25, starting in the top left sub-figure and progressing row-wise to the bottom right sub- 
figure, e.g. t — 2.75 is in the middle column of the top row. 




Time 

Figure 3. The max norm of the solution for < i < 20 for the case b = 0.5, starting from a 
Gaussian pulse. The blue dots correspond to grid size h = 10~^ and the red crosses have 
h = 5x 10-3. 



the surface wave decays slower and slower in the x-direction. In the limit P — 0, the amplitude of the wave is 



28 



TITLE WILL BE SET BY THE PUBLISHER 



constant in x which corresponds to one-dimensional wave propagation in the y-direction, consistent with the 
limiting boundary condition Ux = 0. There are no numerical difficulties in this limit. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 4. The real part of the initial data for the surface wave with j3 — 0.5 and loq ~ Stt. 

The case |/3| — ^ 1 is more difficult to solve numerically. Here we study 0.5 < /? < 1 and we use (|124[) as an 
approximation of the exact solution {u^ is only exponentially small at a; = 1 and does not exactly satisfy the 
boundary condition at that boundary). To make sure the amplitude of the surface wave is negligible at the 
X — 1 boundary, we choose 

ujQ = Stt, e-l'3'^«l = e'^'' « 3.48 x 10"^ /3 = 0.5. 

In Table Owe show the max norm of the error Ug — v for different values of /3. The case (3 = 0.5 shows second 
order convergence, both at time t = 1 and t — 10. As can be expected in wave propagation problems, the error 
is dominated by the phase error, which explains why it is about 10 times larger at i = 10 compared to t = 1. 
For /? = 0.9, the error still converges to second order accuarcy at time t — 1, but shows an unexpected pattern 
at time t = 10. Here the error is larger for the intermediate grid size h — 5 x 10^^ than for the coarse grid size 
h = 10~^. This behavior is explained by studying the time history of the error, see Figure [5l For h = 10~^, 
the max error occurs at t « 5.5 when the numerical solution is about 180 degrees out of phase with the exact 
solution. At later times the error in the numerical solution decreases because it is between 180 and 360 degrees 
out of phase. The grid with h — 5 x 10~^ is barely fine enough to capture the solution at time t = 10 because 
the phase error exceeds 90 degrees. As a result we don't see the expected second order convergence when the 
grid is refined to h = 2.5 x 10"'^. However, the error at t = 10 is about 10 times larger than a.t t ~ 1 for the 
finest grid, which indicates that this resolution is adequate for /3 = 0.9. The situation is even more dire for 
/3 = 0.99. Here the errors at time t = 1 show a simular behavior as at t = 10 for (3 = 0.9, so only the finest 
grid provides adequate resolution at t = 1. At i = 10, the error displays a completely erratic behavior with the 
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largest error for the finest grid. An even finer grid would be necessary to obtain an accurate solution at t = 10, 
when /3 = 0.99. 



Case 


h 


\\Us - v\\co{t = 1) 


\\Us - v\\oc{t = 10) 


P = 0.5 


1 X IQ-'^ 


2.44 X 10"^ 


2.38 X IQ-i 




5 X 10"^ 


6.35 X 10-3 


6.19 X 10-2 




2.5 X 10-3 


1.60 X 10-3 


1.56 X 10-2 


/3 = 0.9 


1 X 10-^ 


6.04 X 10-^ 


2.46 X 10-1 




5 X IQ-^ 


1.58 X 10-1 


1.40 X 10° 




2.5 X IQ-^ 


4.00 X 10-2 


3.95 X 10-1 


/3 = 0.99 


1 X 10"^ 


1.67 X 10" 


1.37 X 10" 




5 X IQ-^ 


5.20 X 10-1 


1.81 X 10-1 




2.5 X 10-3 


1.44 X 10-1 


1.47 X 10° 



Table 2. Max error in the solution as function of the grid size when the exact solution is the 
surface wave Us{x, y, t). 




01 23456789 10 



Figure 5. The max norm of the error as function of time for a surface wave with /3 = 0.9 
computed on a grid with h = 10-2 (blue), h = 5 x 10-3 (green), and h = 2.5 x 10-3 (red). 

So why is it so hard to calculate an accurate numerical solution as \(3\ — > 1? The spatial resolution in terms 
of grid points per wave length only depends on wq. With ujq = 8tt, the wave length is 1/4 and grid sizes 
h = 10"2^ 5 X 10"3j 2.5 X 10"3 correspond to 25, 50, and 100 grid points per wave length, respectively. The 
exponential decay in the x-direction only depends weakly on (3 and never exceeds e"!'^"'^ for (3 < 1. Hence the 
solution varies on the same length scale in the x- and y-directions. Furthermore, the temporal resolution in 
terms of time steps per period only improves as — >■ 1 because the wave speed goes to zero in this limit. We 
conclude that the numerical difficulties are not due to poor resolution of the solution. 

To further analyze the cause of the poor accuracy in the numerical solution for \/3\ — )■ 1, we decompose the 
problem (I108p - (|112p into two parts, 

y,t) = U {x, y, t) + u{x, y, t), 
such that U satisfies a doubly periodic problem on an extended domain, 

Utt = U^^ + Uyy+F{x,y,t), -l<x<2, 0<y<l, t > 0, 
subject to initial conditions. 



C/(x,y,0) = /i(x,2/), C/t(x,y,0) = /2(x,y), -1 < x < 2, < y < 1, 
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and periodic boundary conditions 

U{x,y,t) ^U{x,y + l,t), -1 < a; < 2, t > 0, 
U{x, y, t) = U{x + 3,y,t), < y < 1, t > 0. 

The interior forcing function and the initial data can be smoothly extended to become 3-periodic in the x- 
direction, without changing them on the original domain, 

F{x,y,t) = F{x,y,t), fk{x,y) ^ fk{x,y), < x < 1, < y < 1, t > 0. 

The problem for U is independent of the 6-coefRcient in the boundary condition and can easily be solved 
numerically. 

The difference u' = u — U satisfies the scalar wave equation (I108p - (|112p with homogeneous interior forcing, 
homogeneous initial data, but inhomogeneuos boundary conditions, 

u',~bu'y^g'a{y,t) x = 0, 0<y <1, t>0, (125) 
u',^g[{y,t) x = l, 0<y<l, t>0, (126) 

The boundary forcing functions depend on U according to 

9oiy,t) = 9o{y,t) - {U^{0,y,t) -bUy{0,y,t)) , < y < 1, i > 0, 
g[{y,t)^gi{y,t)-U^{l,y,t), < y < 1, i > 0. 



The corresponding half-plane problems were analyzed in section 12.21 The accuracy problems are unlikely to 
arise from the Neumann boundary condition at x = 1 since it is independent of the 6-coefficient. However, the 
half-plane problem subject to (|125|) satisfies the estimates of Theorem 12.81 Here, b = i /3 corresponds to case 
2), and estimate l\'S9\i shows that the Laplace-Fourier transform of u' satisfies 

|i*'(0,a;,s)|2<-^Mf , Res = r;>0, (127) 
1 — 77^ 



for (w,s) in the vicinity of the generalized eigenvalue sq = ±j -y/l — /S^ wq. In general, the solution becomes 
unbounded as |/3| — > 1. The truncation error terms which perturb the numerical solution are therefore amplified 
by a factor 1/a/1 — /3^, which explains the poor accuracy in the numerical solution as |/3| — 1. 
For boundary data goiv^t) which have a Laplace-Fourier transform that can be written as 

5o(w,s) = sG{u,s), 

estimate (|127p becomes 



\n{0,u.,s)\ <i_^2 ^2 ^|,o|2/^2 ^2 =C(3u, — , s-,so. 



Hence the — ?> 1 singularity cancels out and the solution is bounded independently of (3. The factor 's' on 
the Laplace transform side corresponds to a time-derivative on the un-transformed side. Hence, the solution 
is bounded independently of /3 if the boundary forcing can be written as a time-derivative of a function with 
bounded Laplace-Fourier transform. 



goiy,t)^Gtiy,t), G(y,0) = 0, 



00 



e--^^^^ye~'*G{y,t)dtdy 

2/=0 Jt=0 



< 00, Re s > 0. 
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The latter condition is satisfied if G{y,t) is in L^, i.e., 

1 poo 

/ \G {y, t)\ dtdy <oo. (128) 

y=QJt=0 

To test this theory numericahy, we use a homogeneous interior forcing {F = 0) and homogeneous initial 
conditions (/i = /2 = 0), homogeneous forcing on the x — 1 boundary {gi = 0), and consider three different 
forcing functions on the a; = boundary: gQ^\y,t) = G{y,t), gl^\y,t) — Gt{y,t), and g^\y,t) = Gtt{y,t)- 
Here we choose G{y, t) to trigger a surface wave: 

G(y,t) =u,(0,y,i)e"^*/*''-'^', = 0.2, 

where Ug is defined by (|124p . The Gaussian pulse exp(— (t/to — 7)'^) decays exponentially fast away from its 
center at i = 7io- For example, it equals 1.23 x 10""* at t ^ 7to ± 3to, and 5.24 x 10"^^ at i = 7to ± 7to. The 
function G{y,t) satisfies (|128l) . so our theory predicts that boundary forcings g^^^ and g'"^^ should give solutions 
that are bounded independently of /3. However, the time- integral of a Gaussian pulse is the error- function (erf), 
so the boundary forcing ^q"'^' does not satisfy (|128l) . 

In the numerical calculations we take = Stt and study the cases (5 — 0.5, /3 = 0.9, and /? — 0.99. The grid 
size and time step are h = 2.5 x 10^^ and 5t — 0.5h. The max norm of the solution as function of time is shown 
in Figure [HI The case g^^^ — G in the top sub-figure illustrates the general case where the solution grows as 
1/3 1 1. Note that estimate (|127l) predicts the solution to grow like — which means that the solution 

should be about 3.5 times larger for /? = 0.99 than /3 = 0.9. In the numerical calculation, the max norm of the 
solution grows from about 0.75 for /? = 0.9 to 3.75 for P = 0.99, which is slightly faster than predicted by theory. 

(3) 

The case = Gu in the bottom sub-figure shows the opposite situation when the solution decays as /3 — ?> 1 
because the forcing function is a second time-derivative of a function with bounded norm, corresponding to 
an factor on the Laplace transform side. The intermediate case ^g^'' = Gt is shown in the middle sub-figure. 
Here the solution grows between /3 = 0.9 and (3 = 0.99, but not as fast as for gj^^K To more closely study the 
behavior near /3 = 1, we take /3 = 0.995, 0.999 and 0.9997 corresponding to y/l- ~ 0.0998, 0.0447 and 
0.0244, respectively. To properly resolve the solution we here use an extra fine grid with h = 1.25 x 10~^ and 
St ~ 0.5ft.. The max norm of the solutions, shown in Figure [71 reveal that the solution indeed stays bounded 

(2) 

independently of /3, confirming our theory also for boundary forcing g^ — Gt- 



Appendix 

In this appendix we collect a number of auxilary lemmas. 
Lemma A.l. The solution of 

y^^Xy + F, Re A > 0, < a; < oo (129) 

satisfies the estimate 

Proof. Integration by parts gives us 

iy,y.) ^ -\ym^ - iy..y). i-e., 2Re(2/,y,) = -|2/(0)p. 
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150 



100 




Figure 6. The max norm of the solution as function of time for the boundary forcing functions 
g^^^ = G (top), (7q^^ = Gt (middle), and g'^^ = Gtt (bottom). In each figure, the green, blue, 
and red curves correspond to (3 = 0.5, /3 = 0.9, and (3 = 0.99, respectively. 



Therefore 



-|y(0)p+ReA|lyf = -Re(y,F)<|l2;|l 

With a = 2 the first inequality follows. With a = 1 the second inequality follows. 
Lemma A. 2. The solution of 

yx = ->^y + F, y{0)^g, Re A > 0, < a; < oo, 

satisfies 



□ 



(130) 
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Proof. 



{y,y)x = 2Re {y,y^ 



-2(ReA)|y|2 + 2Re(y,F). 



As ?/ € L^, integrating we have 



-|y(0)|2 = -2ReA \\yf + 2Re{y,F) 
< -2RcA ||y||2 + 2 ||y|| ||F|| 

<-2ReA|lj/f + ReA|b||2 + ^. 



Thus, 



and the lemma follows. 



ReA < |y(0)|2 



ReA 



□ 



Lemma A. 3. Let a,b be real and consider \/a + ib with —it < arg(a + ib) < tt, argy^ = ^argz. Then, the 
following inequalities hold 



2"'/V|a| + IH < I ATl6| < + 
2-^^^\\/a + ib\ < 2"^/''V|a| + \b\ < RcVa + ib < ^a + ib] for a > 0, 
0-5/4 1^1 < 



|%/a + i6| 



^/H + \b\ 



< ReVa + ib < \Va + ib\ for a < 0. 



(131) 
(132) 

(133) 



Proof. In polar notation a + i6 = pe'^, p = Va^ + > 0, — tt < < tt, and 



\/a + ib = \/pe* 



We have 



\/|a| + |6| = VPv'Icosfll + lsin^l > 2i/4y;5 
and the first inequality in (|131l) follows. The second inequality in (|13ip follows from the triangle inequality. If 
a > then | € [-7r/4,7r/4] and (|13ip implies 

l^/^m^l > Re V^TIb = ^cos{ei2) > > 2"^^^^/\a\ + \b\ > 2-^/^\V^^Tib\ 
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which is (|132p . To prove (|133p notice that, as a < 0, 
\b\ ^ \b\ _ psiiK 



V|a| + \b\ ~ \V^Tib\ VP 

therefore 



< , ' ' , = !—— = sin(6l/2)cos(6l/2) < 2 Re Va + ib 

' +ib\ 

1 |6| 



9 ^ <ReVa + ib<Wa + ib\ 

and the inequahty fohows from (|13ip . □ 



For the forthcoming lemmata, we remind the reader that s = r/ + i^, where 77, ^ € R, and k — + s^. We 
shall now apply the last lemma to 



K = y^oj^ + 772 - ^2 ^ 2i£,T], i.e., a = uj'^ + nf - , b = 2(,r]. 

In the following three lemmas we denote by (5 a constant with < (5 < 1. 
Lemma A. 4. Let 

Si^2-^/^VS, <52 =2^/2(1 1/4^ ,53 =min(5i,(52). 
Then 

> { I I (134) 
I i52\/ v I^P + ?7 otherwise. 

Also, always 

\k\ > (5377. (135) 

Proof. By (|13ip we obtain, for the first case, 

\k\ > 2-i/Vl^' + - ^'1 + 2|^|77 > + |s|2. 

If |cj2 + ^2 _ ^2| < J(^2 |^|2)^ ^l^gj^ 

2^^ > (1 - <5)(u;2 + ^2 ^ ^^2) ^ (1 _ ^)(^2 _^ |^|2)_ 

Therefore 



1^1 > 2-1/4^2^ > ^2VT^^A^^TWv- 



Also \/oj^ + |s|2 > 77 implies p35p . This proves the lemma. □ 
Lemma A. 5. 

'2-5/4|k| if UJ^+T]^ -^'^> 0, 



Re K > 

ReK>(54?7, ^4 = 2^^/*min(l,(53). 
Proo/. If a;2 + 7^2 _ ^2 > ^^len p32t gives us 

ReK> 2-3/4|k| > 2-3/%,y. 
If a;2 + 77^ - < 0, then 2^2 > ^^2 ^2 _^ ^2 ^ ^2 _^ |^|2 Therefore ([TSS]) gives us 



Re. > 2-^/4^ > 2-3/4 + 2-3/4^ 
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This proves the lemma. □ 
Lemma A. 6. 

\k\ ReK> Se^yuj^ + \s\^7i 5q = mm(5i(54, 2-"^ 51,2-'^/^'^ . 
Proof. If a;2 + - ^2 > |^2 + ,y2 _ > ^(-^2 ^ |g|2)^ ^-^en and give us 

|k| ReK > (5iv/w2 + |s|2 (54?7. 
If a;2 + 772 - ^2 > g^j^^j 1^2 ^ ^2 _ ^2| ^ j(^2 ^ |g|2)^ i-j^gj^ and ^M)) give us 

|k| ReK > 2-5/'*|k|2 > 2-5/4j2^^2 ^ j^p^. 
Finally, if a;2 + 7^2 - ^2 ^ ([736]) we obtain 

|K|ReK>2-3/V|w|2 + |s|2,7. 

This proves the lemma. □ 

Lemma A. 7. j4ssu7ne t/iat, for the boundary condition 1), a > 0, \b\ < 1. T/ien there is a constant d > such 
that, for all lj and s with Re s > 0, 

\s + aK- ibuj\ > ^^|s|2 + |w|2. 

For the proof, see Lemma 3 of [4]. 

Finally we have a lemma similar to Lemma A. 4 and Lemma A. 6 but for the normalized variables 



k' = y/r]'^ ~ +Lu'^ + 2i^'Lo', C'^ + = 1, |r/| < 1. 

By Lemma A. 3, 



W\ > 2-i/vi-e"+^" + '?"l + 2ie'ih'i 

Re k' > 2-3/4^1 - ^'2 + ^'2 + ^/2| + 2\e\ \r]'\ if ^'^ < u'^ + r]'^, 

Re,c'>^_ l^'l^' if ^'2 >c^'2+,y'2. 

- ^|_^'2+^,2+^,2|+2|e|h'| 

Lemma A. 8. There is a constants S > such that 

\ReK'\>6r]\ \k'\>2-^^%\ \k'\ \Re k'\ > Sev' ■ 
where Sq is that of Lemma A. 6. 

Proof. From Lemma Lemma A. 6, dividing by cj2 + ^2^ 



(136) 



I 'I iRo 'I 1^1 + V . . , 

K ReK = — > de — , — , > Oq ri . 



Now if a;'2 - ^'2 > 0, then w'2 - ^'2 + r]'^ > 7/'2, and 



so that, by Lemma A. 3, 



'|>2-i/V, and RcAi'>2-3/V. 
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If - < then > 1/2. From (fT^ 

>2-i/y2\/2?7'>2i/2 

This last estmiate also holds for |Re«;'| when w'^ — ^'^ + 77'^ > 0. When w'^ — ^'^ + 77'^ < we just use 
|^'2 _ ^'2 ^/2| _^ 2|^'| rj' < const, and from ([TM)) the esthnate for Rere' foUows from ([TM)) . This proves the 
lemma. □ 



We thank Prof. Oscar Reula who participated in the original discussions leading to this work. 
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